Cosmic Rulers 
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We derive general covariant expressions for the six independent observable modes of distortion of 
ideal standard rulers in a perturbed Friedmann-Robertson- Walker spacetime. Our expressions are 
gauge-invariant and valid on the full sky. These six modes are most naturally classified in terms of 
their rotational properties on the sphere, yielding two scalars, two vector (spin-1), and two tensor 
(spin-2) components. One scalar corresponds to the magnification, while the spin-2 components 
correspond to the shear. The vector components allow for a polar/axial decomposition analogous 
to the B/B-decomposition for the shear. Scalar modes do not contribute to the axial (_B-)vector, 
opening a new avenue to probing tensor modes. Our results apply, but are not limited to, the 
distortion of correlation functions (of the CMB, 21cm emission, or galaxies) as well as to weak 
lensing shear and magnification, all of which can be seen as methods relying on "standard rulers" . 



I. INTRODUCTION 

One of the primary goals of cosmology is to accurately 
measure the expansion history and the growth of struc- 
ture in the Universe. Many of the cosmological probes 
used for this purpose can be classified as standard can- 
dles or standard rulers. The most obvious examples 
are type la Supernovae and the baryon acoustic oscilla- 
tion feature in galaxy correlation functions, which in an 
unperturbed Friedmann-Robertson- Walker (FRW) Uni- 
verse directly measure the geometry and expansion his- 
tory of the Universe [iHl]- Beyond the background cos- 
mology, cosmological perturbations affect the apparent 
scale of rulers, which can be used as a probe of structure 
in the Universe. In fact, from this point of view, standard 
rulers comprise a much larger set of observations: for ex- 
ample, galaxy redshift surveys measure the correlation 
function of galaxies, which is then compared with pre- 
dictions based on a cosmological model; in other words, 
the correlation length of galaxies (or any characteristic 
scale in their correlation function) serves as a standard 
ruler. Weak lensing shear, measured using galaxy ellip- 
ticities, uses the fact that galaxies sizes measured along 
fixed directions are on average equal. On the other hand, 
lensing magnification measurements rely on the fact that 
galaxies have a characteristic luminosity (standard can- 
dle) and/or size (standard ruler). Of course, in the latter 
three cases the "ruler" has a large amount of scatter, 
so that one might call it a "statistical ruler". Another 
example of this kind is lensing reconstruction on diffuse 
backgrounds such as the cosmic microwave background 
(CMB) or 21cm emission from the dark ages 0-0] . In this 
approach one uses the intrinsic correlation pattern of the 
background, which is known statistically, to reconstruct 
the distortion from the observed pattern. 



There is a simple, unified description of these various 
cosmological probes: we observe photons from two differ- 
ent directions and redshifts, which correspond to a known 
physical scale (e.g., the comoving sound horizon at re- 
combination, or the characteristic size of a galaxy). In 
this paper, we study in a general covariant setting which 
underlying properties of the spacetime can be measured 
with an ideal standard ruler, working to linear order in 
perturbations. Since we have six parameters to vary 
when scanning over photon arrival directions and red- 
shifts, we can measure six degrees of freedom. These 
can be interpreted as the components of a metric (of 
Euclidean signature) mapping apparent coordinate dis- 
tances into actual physical separations at the source. It 
is useful to further decoixipose these components into 
parts parallel to the line of sight (longitudinal), trans- 
verse, and mixed longitudinal-transverse parts. This is 
equivalent to a decomposition into scalars, vectors, and 
tensors in the two-dimensional subspace perpendicular to 
the photon 4-momentum and the observer's four-velocity; 
i.e. we classify components in terms of their transforma- 
tion properties under a rotation around the line of sight. 
Note that this is independent of the usual decomposi- 
tion of metric perturbations on three-dimensional spa- 
tial hypersurfaces (i.e., in terms of the transformation of 
plane- wave metric perturbations under a rotation around 
the fc-vcctor). We will denote the latter ("3-scalars" and 
so on) as C, Bi, Aij, and the former ("2-scalars" etc) 
as C, Bi, A' 
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The transverse components of the dis- 
tortion, Aij, are perhaps best known. They correspond 
to the magnification (2-scalar) and shear (2-tensor). We 
show that in general one can also measure a longitudinal 
scalar and the two components of a vector on the sphere. 

On scales much smaller than the horizon, an effective 
Newtonian description is sufficient, and this is what es- 
sentially all previous studies are based on. However, up- 
coming surveys will probe scales approaching the hori- 
zon, and an iirterpretation of these data sets can in prin- 
ciple be hampered by gauge ambiguities. In the case 
of the correlation of galaxy density contrast, this issue 
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has attracted significant interest and has recently been 
resolved [8l-[TT|. The unified treatment presented here 
resolves these issues for the wide set of cosmological ob- 
servables mentioned above. More precisely, we obtain 
general coordinate-independent and gauge-invariant re- 
sults for all observables, including the shear and mag- 
nification. The differential equation {optical equation) 
governing the magnification and shear was first derived 
in [13]. The magnification has been derived to first or- 
der in [ll]. The shear has been derived to second order 
in conformal-Newtonian gauge in [ij], while [TH| derive 
the shear for general backgrounds. To the best of our 
knowledge, the expression for the observable shear writ- 
ten in a general gauge is presented here for the first time. 
Further, all expressions are valid on the full sky. Our 
approach naturally includes the "metric shear" contribu- 
tion introduced in [l6| , and we provide a straightforward 
physical interpretation of our result. 

In addition, we show how the (2-)vector observable un- 
covered here can be decomposed into E- and i3-modes in 
analogy with the shear, corresponding to polar and axial 
vector parts. As in the case of shear and CMB polar- 
ization, (3-)scalar perturbations do not contribute to the 
_B-mode, while (3-)tensor perturbations contribute. This 
in principle offers another avenue to search for a stochas- 
tic gravitational wave background in large-scale struc- 
ture, since no scalar perturbations contribute at linear 
order. However, a spectroscopic data set is likely neces- 
sary to reconstruct the vector component with an inter- 
esting signal-to-noise. 

Apart from the linear treatment of metric perturba- 
tions, we make two further simplifying assumptions: first, 
we assume "small rulers" in the sense that rulers subtend 
a small apparent angle and redshift interval. Wide-angle 
effects are likely negligible for almost all applications (the 
large-scale BAO feature being perhaps the most impor- 
tant exception) . A treatment of wide-angle effects neces- 
sarily involves a detailed model of the survey geometry, 
which is clearly beyond the scope of this paper. The sec- 
ond assumption is that any scatter or variation in the 
actual ("intrinsic") physical scale of the standard ruler 
is uncorrelated with large-scale perturbations. This will 
not hold true in general, since the physical systems used 
as rulers will be affected by their large-scale environment. 
One well-known example is the intrinsic alignment con- 
tribution to shear correlations [l7| . Further examples 
include the distortion of correlation functions by large- 
scale tidal fields [3, or by a non-Gaussian coupling of 
the density field to primordial degrees of freedom [l9[. 
Since these "intrinsic effects" depend on the physics of 
the given ruler, we refrain from discussing them here, as 
they would distract from the generality of the rest of the 
results. 

Finally, while we focus on general standard rulers here, 
the case of standard candles is directly related to our re- 
sults. This is because the relation between angular di- 
ameter distance Da and luminosity distance D^, 



where z is the observed redshift, holds in a general space- 
time and for any source (this is a consequence of pho- 
ton phasespace conservation). Thus, the magnification 
measured for standard candles is identical to the mag- 
nification for standard rulers which we will derive here. 
However, as discussed, standard rulers can measure five 
additional degrees of freedom not accessible to standard 
candles. 

Our results are of immediate relevance to recent stud- 
ies which consider the i?-modes of the cosmic shear as 
possible probe of an infiationary gravitational wave back- 
ground [ll,[20,[2l|, and to studies that propose to use the 
high-redshift 21 cm emission for the same purpose (22l . [23j . 
In particular, our expressions can be used directly to con- 
struct optimal estimators on the full sky for searching for 
the imprint of gravitational waves in a three-dimensional 
field such as the 21cm background. We also use many of 
the results derived here in two recent papers studying the 
impact of gravitational waves on the observed large-scale 
structure [Si, HI. 

The outline of the paper is as follows: we begin in § HI] 
by introducing our metric convention and useful nota- 
tion. The general expression for the mapping from ap- 
parent size to true physical size of the ruler is derived in 
§ mil We then decompose the contributions into longitu- 
dinal and transverse parts in § IIVI The following three 
sections deal with these different parts consecutively. We 
discuss and conclude in § IVIIII The appendix contains a 
large amount of additional reference material on multi- 
pole expansions of higher spin functions, perturbed pho- 
ton geodesic equation, and various test cases applied to 
our results. 



II. NOTATION 



by 



In a general gauge, the perturbed FRW metric is given 



ds^^ 0^(77) - {I + 2A)dr]^ - 2B,d7jdx' 



+ (Sij + hij) dx^dx-' 



(2) 



where we have assumed a spatially flat Universe (curva- 
ture can be included straightforwardly, at the expense of 
some extra notation). Here, 77 denotes conformal time. 
Often, the spatial part is further expanded as 



2D5i. 



2E,. 



(3) 



Dr 



{1 + S)'Da 



(1) 



where Eij is traceless. We shall also present the most in- 
teresting results in two popular gauges: the synchronous- 
comoving (sc) gauge, where A = = Bi, so that 

ds"^ 0^(77) [-dif + {Sij + hij) dx'dx^] ; (4) 

and the conformal-Newtonian (cN) gauge, where Bi = 
= Eij. In the latter case, we denote A = 5*, Z) = $, 
conforming with standard notation, so that 

ds^ = a^{r]) [-(1 + 2^)d7f + (1 + 2^)6ijdx'dx^ . (5) 
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We also denote the background FRW metric (in the ab- 
sence of perturbations) as g^^ = o?{ri)ri^^. 

It is useful to define projection operators parallel and 
perpendicular to the observed line-of-sight direction fi*, 
so that for any spatial vector X"^ and tensor Eij , 

X\\ =n,X\ 

V'^=5'^-fenK (6) 

Correspondingly, we define projected derivative opera- 
tors, 

9|| = n*9i, and 

d'^=V''dj. (7) 

Note that 9^, 9|| and d\, d''^ do not commute. Further, 
we find 

where x is the norm of the position vector so that = 
Ix- Note that and dn commute. More expressions 
can be found in § II of 

Finally, it proves useful to decompose the quantities 
defined on the sphere, i.e. as function of the unit line- 
of-sight vector n, in terms of their properties under a ro- 
tation around fi. In particular, consider an orthonormal 
coordinate system (ei, 62, n). If we rotate the coordinate 
system around n by an angle ?/', so that — >• e-, then 
the linear combinations m± = (ei =p ie2)/-\/2 transform 
as 

m± ^m'± ^e^^'^mi. (9) 

We say that a general function /(fi) is spin-s if it trans- 
forms under the same transformation as 

/(n) ^ /(n)' = e^-^/(n). (10) 

An ordinary scalar function on the sphere is clearly spin- 
0, while the unit vectors m-i- defined above are spinil 
fields. More details can be found in App. [Xj This de- 
composition is particularly useful for deriving multipole 
coefficients and angular power spectra. We also define 

E± = m'^m^^E.j (11) 

for any 3- vector Xi and 3-tensor Eij . 

For the quantitative results shown in Fig.[3l we assume 
a flat ACDM cosmology with h = 0.72, 17„ 0.28, a 
scalar spectral index Us ~ 0.958 and power spectrum 
normalization at z = of cts = 0.8. 



III. STANDARD RULER 

In the absence of perturbations, photon geodesies are 
given by straight lines in conformal coordinates, 

x^'{x)^{Vo-X,nx), (12) 

where we have chosen the comoving distance x affine 
parameter. Correspondingly, for a photon arriving from 
a direction n with redshift z, we assign an "observed" 
position of emission x'^ given by 

i° = 770 - X 
=n'x 

x = m, (13) 

where x(5) denotes the comoving distance-redshift rela- 
tion in the background Universe. The actual spacetime 
point of emission, denoted with x^^ , is displaced from the 
observed positions by Ax^ (see also Fig. [T]), 

x'^ = 5^ + Aa;'^(n,i). (14) 

Further, we will need the scale factor at emission. It is 
related to the inferred emission scale factor a = 1/(1 + z) 

by 

a(x'^(h, z)) . , . ^ X , ^ 

— = 1 + Alna(n,f). (15) 
a 

At first order, Alna = Az/(l-|-z), where Az is the differ- 
ence between the observed redshift and the redshift that 
would be observed in an unperturbed Universe. Note 
that the latter quantity is gauge-dependent. We will give 
explicit expressions for A In a and Ax'' in § IIVI 

The displacements Alna, Ax'' are not observable 
(they depend on which gauge, or frame, the spacetime 
perturbations are described in). In order to determine 
actual observables, we consider the case of a standard 
ruler. A standard ruler exists if we can identify two 
spacetime points which are separated by a fixed space- 
like distance rg. What we observe is the apparent size 
at which this ruler appears in a given direction n and 
redshift z. Let fi, z and n',z' denote the observed coor- 
dinates of the "end points" of the ruler, and x and x' 
the apparent spatial positions inferred through Eq. . 
The inferred physical separation is then given by 

f2 ^ a^S,j{i' - S:"){S:^ - i'^), (16) 

where a ~ 1/(1 + z) is the observationally inferred scale 
factor at emission (Fig. [T|). 

We now have to carefully consider what the condition 
of a standard ruler in cosmology means. A useful, physi- 
cally motivated definition is that it corresponds to a fixed 
spatial scale as measured by local observers which are co- 
moving with the cosmic fluid; precisely, the spatial part 
of the four-velocity u'^ of these observers is given by 

V' = (17) 
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a{x°) 
a - + 



Oo = 1 



FIG. 1: Illustration of the apparent and actual standard ruler. 
Photons arrive out of the observed directions n, n' and with 
The apparent positions are indicated 



by X 



^'^ x''^ , while the true positions are at x'' 



X perturbed 



by the displacements Aa;'', Aa;''' (whose magnitude is greatly 
exaggerated here), f is the apparent size of the ruler, while 
To is the true ruler. 



We are mostly interested in applications to the large-scale 
structure during matter domination; in this case, the cos- 
mic fluid is simply matter (dark matter -t- baryons) , and 
there is no ambiguity in this definition. In synchronous- 
comoving gauge, Eq. (|17|) yields — 0. Further, for now 
we assume the ruler is fixed. We consider an evolving 
ruler in § 11101 

This definition can also be phrased as that the length 
of the ruler is defined on a surface of constant proper time 
of comoving observers. This proper time corresponds to 
the "local age" of the Universe. The separation of the 
two endpoints of the ruler, x*^ , x'^, projected onto this 
hypersurface should thus be equal to the fixed scale tq: 

+ u,,{x")u,{x")] {x^^ - x'^)(x^ - x'^ - r^, 

(18) 

where -I- is the metric projected perpendicular 

to Ufj_, the four- velocity of the comoving observers (note 
that u^u^ = —1). Here and throughout, we will assume 
for simplicity that the ruler is "small", i.e. it subtends 
a small angle, and redshift interval {\z — z'\ z). This 
entails f ^ x, and that wc can simply evaluate the metric 
and four-velocity at either end-point (corrections involve 
higher powers of x^ — x'^). 

The four- velocity of comoving observers, whose spatial 
components are fixed by Eq. (jl7p . is given by 



1 - A, v' 



- a (-1 - A, Vi- Bi) , 



(19) 



where we consider to be first order (as the metric per- 
turbations). In the following, we will assume sources to 
be comoving as well, i.e. to follow Eq. ([T^ . It is straight- 
forward to generalize the treatment to different source 
velocities. Using Eq. ([2]) and Eq. (fT9|) . wc have 



( 



V 



(20) 



With this, Eq. dH]) yields 

- 2~c?v^ {Ji^JF -f 5£°[Aa;' - Ax'*] + (5F[Ax" - Ax-'"]} 
+ g„(x")|(5iMi^' + 6i\l\x^ - Ax'J] + [Ax* - Ax'*],5i^'| 

(21) 

where Ax^ = Ax^(n, z). Ax''' = Ax^(n',z'), and the 
components of the apparent separation vector are 



(22) 



In order to evaluate the spatial metric gij{x°') at the 
location of the ruler, we use Eq. ([15]) to obtain at first 
order 



g^^ (x") = [(1 + 2A In a) % + h. 



(23) 



We now again make use of the "small ruler" approxima- 
tion, so that 



Ax* - Ax'* 



d 

5£"— Ax\ 
9x" 



(24) 



Like any vector, we can decompose the spatial part of the 
apparent separation 5x* into parts parallel and transverse 
to the line of sight: 



(5x|| = fiiSx^ 

Sx'^ EE V'jSi^ = Si' - h'Sin 



(25) 



In the correlation function literature, ^X||, |(5xj_| are 
sometimes referred to as ir and a, respectively. Then, 



d 



Sx"— = {Sx^'d,^ + (5.^11 9||) + Sx'^d^, 



(26) 



where wc have similarly defined d\\ — Wdi, dj_i = V^'' dj. 
Since the observed coordinates x'' by definition satisfy 
the light cone condition with respect to the unperturbed 
FRW metric, we have 5x^ ~ "'^^11 small-angle 
approximation. Thus, 

Sx'^dr; + <^5||9|| = (5511(911 — dri) 

^6S:\S=Si\\HCz)-^. (27) 
" ox oz 

where d/dx is the derivative with respect to the affine 
parameter at emission. We thus have 



i9x" 



(5x||9^ -I- 5x\d±i. 



(28) 
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Working to first order in perturbations, we then obtain 



M 



rl-f^^ 2Alna + a^hijSx^Sx^ 

+ 2a^(u||(5x| + v^i5x\5x\\) 

+ 2h^5^j5x' (<5.i ii^x + 5x\dA^ k) Aa;^ (29) 

All terms are straightforward to interpret: there are the 
perturbations to the metric (both from the metric per- 
turbation itself and the perturbation to the scale factor 
at emission); the contribution cx v from the projection 
from fixed-ry to fixed-proper-time hypersurfaces; and the 
difference in the spatial displacements of the endpoints 
of the ruler. 



A. Evolving ruler 

So far, we have assumed that the physical scale of the 
ruler is fixed. This does not have to hold in general; for 
example, the BAO scale is fixed in terms of comoving 
coordinates. We now consider the case where the (mean) 
physical scale rg evolves over cosmic time. It is simplest 
to consider tq as a function of scale factor; one can easily 
convert to other variables such as conformal time using 
the relation with the scale factor in the background Uni- 
verse. 

Then, on the left-hand side of Eq. ((29| is to be 
evaluated for the scale factor at emission: 



rl[a{x')]^rl{h) 



1 



2 — A in a 

a m a 



(30) 



Thus, if we compare the apparent size of the standard 
ruler to the true size at the apparent time of emission 
(of course assuming that we are able to predict ro(a)), 
Eq. ((^ becomes 



~2 

r 



2Alna 



din ro 
9 In a 



a^hijSx^Sx-' 



+ 2a^(w||(5i| + v± iSx''^6x\\ ) 

+ 2d.^d^jSi' ((5i||% + Si'ld^k) A.T^'. (31) 

Note that if ro oc a, i.e. if the ruler corresponds to a fixed 
comoving scale, the terms multiplying A In a cancel. This 
is as expected, since a perturbation to the scale factor at 
emission does not affect a fixed comoving scale. 




FIG. 2: Illustration of the distortion of standard rulers due 
to the longitudinal (2-)scalar C, (2-)vector B, and transverse 
components, magnification A4 and shear 7. The first row 
shows the projection onto the sky plane, while the second 
(third) row show the projection onto the line-of-sight and 
{x\^) axes, respectively. In case of B and 7, we only show one 
of the two components. See also Fig. 3 in 



case for photometric galaxy surveys. On the other hand, 
spectroscopic surveys and redshift-resolved backgrounds 
such as the 21cm emission from high-redshifts are able to 
measure the longitudinal displacements as well. 

Noting that = a^[5a;J + ((5x^)^], and taking the 
square root of Eq. (|3ip . we obtain the relative perturba- 
tion to the physical scale of the ruler as 



{Sill? 



C 



Sx^^Sx^^ 



5x\5x''^ 



(32) 



where we have defined fc = f/a as the apparent comoving 
size of the ruler. The quantities multiplying C, Bi, Aij 
are thus simply geometric factors. The coefficients are 
given by 



C 

B^ 



A In a 



1 



9 In a 



-\H-n 



d^Ax\ 



- Alna 



1 - 



din ro 
9 In a 



h''d^.,Axk - d^Axj_i 
1 



IV. SCALAR-VECTOR-TENSOR 
DECOMPOSITION ON THE SKY 

It is useful to separate the contributions to Eq. ([^ 
in terms of the observed longitudinal and transverse dis- 
placements. For some applications, only the transverse 
displacements are relevant. This is the case for dif- 
fuse backgrounds without redshift resolution, such as the 
CMB or the cosmic infrared background, and largely the 



-iPJk^^^+V^k^^J)Ax'' 



(33) 



where Aa;|| , Ax^j_ are the parallel and perpendicular com- 
ponents of the displacements Ax*. Note that while we 
have assumed that the ruler is small (i.e. 6x'^ <C x), 
the expressions for C, Bi, Aij are valid on the full sky. 
Fig. [5] illustrates the distortions induced by these com- 
ponents. Observationally, we have 6 free parameters (as- 
suming accurate redshifts are available): the location of 
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one point fi, z, and the separation vector described by 
(with 5x^ being fixed by the hght cone condition). 
Using these, we can measure a (2-)scalar on the sphere, 
C, a 2 X 2 symmetric matrix, Aij^ and a 2-component 
vector on the sphere, Bi. As a symmetric matrix on the 
sphere, Aij has a scalar component, given by the trace 
M. = V^^ Aij {magnification), and two components of 
the tracelcss part which transform as spin- 2 fields on the 
sphere {shear ±27 as defined in Eq. ((57|) below). These 
quantities are observable and gauge-invariant, while any 
of the individual contributions in Eq. p3p are not in 
general. Note that we cannot measure any of the anti- 
symmetric components, such as the rotation. This is be- 
cause we have not assumed the existence of any preferred 
directions in the Universe. If there is a primary spin-1 or 
higher spin field, such as the polarization in case of the 
CMB, then a rotation can be measured as it mixes the 
spin±2 components (see, e.g. [1^). In the next sections 
we study these three terms in turn. 

For reference, we now give the explicit expressions for 
the displacements Aa;' and A In a. They are defined such 
that Ax* = = Alna for a local source, i.e. for z k, Q. 
The details of the derivation are presented in App. |B] 
Separating into line-of-sight and transverse parts, we 
have 



Azii 



Axl = 



A- B« 



1 



^Alna (34) 
X (35) 



+ T^{d\_hjk)h^h'' 



l^V^'{h,k)oh^ + Bl,-v\, 



(36) 



dx 



-{B\ 



+ {x~x)dl[A-Bn--h 



The perturbation to the scale factor at emission is given 
by 



A In a = - A + ui] 



dx 



A~Bn 



1. 



(37) 



Here, a subscript o indicates quantities evaluated at the 
observer, while primes denote derivatives with respect to 
T]. Note the appearance of the scalar quantity A — B^^ — 

inEqs. This is the "lensing potential" 

^ in conformal-Newtonian gauge, written in the general 
gauge Eq. 



In particular, in the two popular gauges introduced in 
§ imi Eq. (1271) becomes 

(Alna),e = ^ rdxh\, (38) 
^ Jo 

(Alna)cN = ^-o - * + i^ii - i^iio + f dx - • 

Jo 

(39) 

The latter result clearly shows the "Sachs- Wolfe" , 
"Doppler" , and "integrated Sachs- Wolfe" contributions. 

V. LONGITUDINAL SCALAR 

The longitudinal component can be simplified to be- 
come 



C 



A In, 



91nro 

dz \H{z) J 9 In a 



A-V11+ Bu 



The first line contains the contributions due to the 
fact that the scale factor at emission is perturbed from 
1/(1 -I- z): first, the evolution of the physical standard 
ruler with scale factor, 1 — 91nro/91na, and second, 
the evolution of the distance-redshift relation. The sec- 
ond line contains the perturbations from the metric at the 
source location {—A) and the projection from coordinate- 
time to proper-time hypersurfaces (i?|| — ■^n)- Finally, the 
contributions from the line-of-sight derivative of the line- 
of-sight displacements (oc (1 + z)/H{z)) are given in the 
third line. Note the term which is the dominant 

term on small scales in the conformal-Newtonian gauge. 
This term is also responsible for the leading-order red- 
shift distortions [2^. Apart from the pertubation to the 
scale factor at emission, C does not involve any integral 
terms; this is expected since C is the only term remaining 
if the two lines of sight coincide (n = fi'). In this case, 
the two rays share the same path from the closer of the 
two emission points, and no quantities integrated along 
the line of sight can contribute to the perturbation of the 
ruler. 

Restricting to the synchronous-conioving and 
conformal-Newtonian gauges, respectively, we obtain 



(C) 



cN 



- (Alna)s 

2H{i) II 

- (Alna)cN 



dz \H{z) 



9 In To 



9 In a 

(41) 
9 In tq 



- - V 



H{~z) 



dz \H{z) J 9 In a 
^ (-9,1* + 9||.|| -.[,+$'). 

(42) 
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FIG. 3: Angular power spectra of the different standard ruler 
perturbations produced by a standard scale-invariant power 
spectrum of curvature perturbations: C, i5-mode of Bi, E- 
mode of the shear, and magnification A4. All quantities are 
calculated for a non-evolving ruler and a sharp source redshift 
of 2 = 2. For comparison, the thin dotted line shows the 
angular power spectrum at z = 2 of the matter density field in 
synchronous-comoving gauge. Note that all quantities shown 
here, except for are gauge-invariant and (in principle) 
observable. 



Note that in case of the sc-gauge expression, the redshift- 
space distortion term is included in the last term, through 
/iJ|/2 = D' + Fig. [3] shows the angular power spec- 

trum of C due to standard adiabatic scalar perturbations 
in a ACDM cosmology (the details of the calculation are 
given in App. [E|. Clearly, C is of the same order as the 
matter density contrast in synchronous-comoving gauge 
on all scales. In particular, the velocity gradient term 
dominates over all other contributions. Due to the dif- 
ferent dependence on the angle with the line of sight, 
the projection kernel of C is proportional to d^ji{x), 
while that of is oc ji{x). The former favors larger 
a; at a given I, and thus leads to a relative suppression 
as the slope of the matter power spectrum changes at 
k > 0.01 h/Mpc. 



VI. VECTOR 



Next, we have the two-component vector 



H{S) 



rrS^jAlna, 



(43) 



where we have inserted projection operators for clarity 
(these are trivial since Bi is contracted with Sx\^). As 
expected; this vector involves the transverse derivative of 
the line-of-sight displacement and the line-of-sight deriva- 
tive of the transverse displacement. Note that these two 
quantities are not observable individually. 

Using the spimbl unit vectors m±, Bi can be decom- 
posed into spimbl components: 



B, = +iBmt 



-iBrnl 



±iB = m'^B^ = ~v± + B± + j^d±A\na, (44) 

where we have used the notation of Eq. (jlip . Similar 
to before, we can specialize this general result to the 
synchronous-comoving and conformal- Newtonian gauges: 



(±lS),c 
(±l^)cN 



1 + S 



2H{z 



/ dx-d±h'u 
z) Jo X " 



(45) 



I + z 
v± H -—9-1- A In a 



v± 



H{~z) 

l + z 
H{~z) 



\\o\ 



\- / dx^a±($'-*' 

Jo X 



(46) 



On small scales, the dominant contribution to Bi comes 
from the transverse derivative of the line-of-sight compo- 
nent of the velocity d±v\\, which is of the same order as 
the tidal field. 

Applying the spin- lowering operator B to iB (see 
App. 1^ yields a spin-zero quantity, which can be ex- 
panded in terms of the usual spherical harmonics^. We 
then obtain the multipole coefficients of B as 



{I 



Y^ Jdn [8iS(n,2)]y4(n). 



(47) 



An equivalent result is obtained for In general, 

the multipole coefficients are complex, so that we 
can decompose them into real and imaginary parts. 



BB 



(48) 



^ This is of course equivalent to expanding iB in terms of spin-1 
spherical harmonics. 



8 



One can easily show (App. [X| that under a change of 
parity afj^ transform as the spherical harmonic coeffi- 
cients of a vector (parity-odd) , whereas afj^ , picking up 
an additional minus sign, transforms as those of a pseudo- 
vector (parity-even) . These thus correspond to the polar 
( "E" ) and axial ( "B" ) parts of the vector Bi . 

As required by parity, scalar perturbations do not con- 
tribute to the axial part afj^ (this is shown explicitly in 
App. |D]). Thus, a measurement of the vector component 
Bi of standard ruler distortions offers an additional pos- 
sibility to probe tensor modes with large-scale structure, 
as tensor modes do contribute to afj^ fApp. lPj) . Thus, in 
principle the axial component of Bi could be of similar 
interest for constraining tensor modes as weak lensing B- 
modes [2l|, though one likely requires accurate redshifts 
to measure Bi to sufficient accuracy. We leave a detailed 
investigation of this for future work. 

The power spectrum of the E-mode of B due to stan- 
dard scalar perturbations is shown in Fig. [3] (see App.IE|. 
While the dominant contribution to C is oc fc^/fc^ (^toII^i 
for a given Fourier mode of the matter density con- 
trast in synchronous-comoving gauge (§ IV]), the cor- 
responding contribution to i3 is cx A:j^fc||/fc^ (5^j(k, z). 
Even though approximate scaling arguments suggest that 
Cc{l), should scale roughly equally with /, we see 

that Cb(/) scales faster with / for ^ < 500. The rea- 
son is that the projection kernel for the £'-mode of B 
(oc {dxji)/x) is relatively suppressed with respect to that 
of C (oc d^ji) at large x/l. Since / < 500 corresponds 
to a typical k < 10~^ h/Mpc at the source redshift, 
where Pm{k) oc fc, larger x/l arc favored for progressively 
smaller /, leading to a more rapid decrease of Cb{1) to- 
wards smaller This suppression is thus fundamentally a 
consequence of the shape of the matter power spectrum. 



system (e6/,e0,n), we can thus write 

. _/X/2 + 7i 72 
"I 72 M/2 - 71 



(50) 



Below, we will derive magnification and shear without 
reference to a fixed coordinate system. 



A. Magnification 



Taking the trace of Eq. yields 



2Alna 



d In 7'o 
9 In a 



-\{h\-h\{)+2k-'^^x\\ 
(51) 



The magnification is directly related to the fractional per- 
turbations in distances (see [13, [H]) through 



AD 



_ ^ AA4 ^ 1^ 
Dl Da 2 ' 



(52) 



where the first equality for the luminosity distance follows 
from Eq. ([!]). The contributions to the magnification are 
straightforwardly interpreted as coming from the conver- 
sion of coordinate distance to physical scale at the source 
(from the perturbation to the scale factor A In a and the 
metric at the source projected perpendicular to the line 
of sight, — /i||); from the fact that the entire ruler is 
moved closer or further away by Aa;||; and finally from 
the coordinate convergence k defined through 



k = -i5_LiAa;j 



(53) 



VII. TRANSVERSE TENSOR: SHEAR AND 
MAGNIFICATION 

Finally, we have the purely transverse component. 



A In a 



d In 7'o 



9 In a 
1 



X 



(49) 



where we have again inserted projection operators for 
clarity (note that Vij serves as the identity matrix on the 
sphere). As a symmetric matrix on the sphere, Aij has 
a scalar component, given by the trace A, and two com- 
ponents of the traceless part which transform as spin-2 
fields on the sphere. The trace corresponds to the change 
in area on the sky subtended by two perpendicular stan- 
dard rulers. Thus, it is equal to the magnification Ai (see 
also Fig. [5]). The two components of the traceless part 
correspond to the shear 7. If we choose a fixed coordinate 



This term dominates the other contributions to A4 on 
small scales. However, the coordinate convergence is a 
gau ge- dependent quantity; see for example App. B2 in 
prj . For the general metric Eq. ([2]) it is given by 



1 



((/i',)o-3(/i||)„) -2(i?|| -7;||)„ 



(54) 



dx 



X 

+ i(/7%-3/i||) 

A 

X I "2 



In conformal-Newtonian gauge, it assumes its familiar 
form. 



(k)c 



1 



«l|o + i / dx?(x-x)Vi (*-$), (55) 
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with an additional term —v\\o contributing to the dipole 
of K only, which corresponds to the relativistic beam- 
ing effect at linear order. Explicit expressions for the 
magnification in different gauges are straightforward to 
obtain, however they become lengthy. Here we just show 
that we recover the result obtained in a different way in 
Jcong et al. [Tf| in synchronous-comoving gauge. Using 
Eq. ©, {h\ - h\\)/2 = 2D ~ E\i. This, and defining 
Sz = (Alnaj.c [Eq. ^\ yields 

/ a In rn \ 2 
(X)sc = -21- — — ^ (5z - 21? + Si, + 2k - - Axii . 
\ olna / " X 

(56) 

For a non-evolving standard ruler ((91nro/91na = 0), we 
thus recover the covariant magnification, M. = bM., as 
derived in [IT|. In general, the observable magnification 
depends on the rcdshift evolution of the standard ruler 
considered; for example, galaxy sizes which can be used 
as standard ruler to measure magnification [29j in general 
show a non-trivial rcdshift evolution. 



B. Shear 

We now consider the traceless part of Aij , given by 

7y (n) EE - \^^JM 

- 9^(jAxj_j) - Vijk. (57) 

Here, the terms oc Vij in Eq. (|49p drop out. The terms in 
the second line here is what commonly is regarded as the 
shear, i.e. the trace-free part of the transverse derivatives 
of the transverse displacements. The first term on the 
other hand is important to ensure a gauge-invariant re- 
sult. This is the term referred to as "metric shear" in [l6j . 
Its physical significance becomes clear when constructing 
the Fermi normal coordinates for the region containing 
the standard ruler. 

Consider a region of spatial extent i?, say centered on 
a given galaxy, with R assumed to be much larger than 
the scale of individual galaxies. We can construct or- 
thornomal Fermi normal coordinates [sO, [Ml around the 
center of this region, which follows a timelike geodesic, by 
choosing the origin to be located at the center of the re- 
gion at all times, and the time coordinate to be the proper 
time of this geodesic. The spacetime in these Fermi co- 
ordinates {tp, x'p) then becomes Minkowski, with correc- 
tions going as x'^/R'^ where Rc is the curvature scale of 
the spacetime. Thus, as long as these corrections to the 
metric are negligible, there is no preferred direction in 
this frame, and the size of the standard ruler has to be 
(statistically) independent of the orientation. The most 
obvious example is galaxy shapes, which are used for cos- 
mic shear measurements. In the Fermi frame, galaxy ori- 
entations are random. Note that the Fermi coordinates 



are uniquely determined up to three Euler angles. The 
statement that galaxy orientations are random in this 
frame is thus coordinate-invariant. 

As an example, consider the case where we have a 
purely spatial metric perturbation (cf. Eq. ([2])) at a fixed 
time. We can then expand around the origin, 

/iy(x) = %(0) + /iy,fc(0)x'=. (58) 

Higher order terms are suppressed by (x/Rc)'^. Now, 
consider coordinates given by 

a-^x'p = x' + lh.M^^ + \l2h,^,^)-h,,M^^x\ 

(59) 

In these coordinates, the metric becomes 

= Vt.. + Oixl). (60) 

Thus, it is in terms of the coordinates x'p that galaxies 
should be isotropically oriented on average, not in terms 
of the cosmological coordinates x' . Correspondingly, in 
order to obtain the shear relative to the Fermi frame, we 
need to add the transformation Eq. ([59|) to the displace- 
ments Aa;*: 

A.T* ^ Ax' + ^h,jiO)x^ + i [2/i,j- fe(0) - hjkAO)] x'x''- 

(61) 

With these new displacements, the transverse derivative 
of the transverse displacement becomes 

(62) 

where the last term is suppressed by the size of the ruler 
over the wavelength of the metric perturbation, and is 
thus negligible in the small-ruler approximation. We 
see that Eq. (p^ agrees exactly with the result derived 
above, Eq. ([57| (after subtracting the trace of Eq. (p^ ). 
In other words, the shear derived in the standard ruler 
formalism (§ IIIip is equivalent to the statement that the 
ruler is isotropic in its Fermi frame, the additional term 
coming from the transformation from global coordinates 
to the local Fermi coordinates. This additional term was 
introduced in [l6j as "metric shear" , with a similar mo- 
tivation as given here. In our case, this term is naturally 
included in the standard ruler formalism. 

^ij is a symmetric trace-free tensor on the sphere, 
and can thus be decomposed into spin±2 components 
(in analogy to the polarization of the CMB). Following 
App.El(see also [S^) we can write 7ij as 

lij = 27 m\m\ + _27 m!_mi 
±21 = ml-rn-l^lio , (63) 

where ±27 are spin±2 functions on the sphere (in analogy 
to the combination of Stokes parameters Q ± iU). We 
obtain for the shear components 
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±27 = - T;f^± - Ax^j 



11 

^h±- - {h±)o- dx 



1-24 [m'^d±Bk + {d±hik)m':^n''] - -h± 



(64) 



Eq. (|64p is valid in any gauge. We can now specialize to the synchronous-conioving (sc) and conformal-Newtonian 
(cN) gauges: 



[±2 



7)s 



11 n 



dx 



1-2^] {d±hM)m%n' - yh± 



1 

J 

X 



(65) 



+ {X- X)^^{m%m^d^djhik)h'h'' 



(±27) 



cN 



dx(x-x)^«a,9, (*-$), 



(66) 



In case of the cN gauge, we have used that hij = 2^Sij, 
and thus h± = 0. We see that Eq. ([M]) recovers the 
"standard" result; in other words, there are no addi- 
tional relativistic corrections to the shear in this gauge. 
This is not surprising following our arguments above: in 
conformal-Newtonian gauge, the transformation Eq. (|59p 
from global coordinates to the local Fermi frame is 
isotropic since hij = 2^Sij. Thus, it does not contribute 
to the shear. Note however that only scalar perturba- 
tions are included in this gauge; when considering vec- 
tor or tensor perturbations, one has to use a different 
gauge, for example synchronous- comoving gauge (see [2l| 
for a study of tensor perturbations). Thus. Eq. (j64p and 
Eq. (|65p are important new results. 

In App. \C\ we apply several test cases to the shear in 
synchronous-comoving gauge, Eq. (j65p . in order to ver- 
ify that it is gauge-invariant and correctly reproduces 
known results. In particular, we consider a Bianchi I 
cosmology which induces a shear due to the anisotropic 
angular diameter distance. We also show that Eq. 
when restricted to scalar perturbations, does not produce 
B-mode shear. 

Fig. [3] shows the angular power spectrum of shear and 
magnification due to scalar perturbations for a sharp 
source redshift z = 2 (see App.|E|. For I > 10, the results 
follow the familiar relation Cm{1) = ^C^^{1), valid when 
all relativistic corrections to the magnification become 
irrelevant so that ~ 2k. These corrections slightly 
increase the magnification for small We also sec that 
7 and M. are suppressed with respect to C and B (on 
smaller scales) , at least when the latter are evaluated for 
a sharp source redshift. This is a well-known consequence 
of the projection with the broad lensing kernel, leading 
to a cancelation of modes that are not purely transverse 



(see e.g. [33|). 



VIII. DISCUSSION 

Over the past decade, cosmology has benefited from a 
vast increase in the available data, which have been ex- 
ploited through a broad variety of methods to constrain 
the history of structure in the Universe. Clearly, this calls 
for a rigorous investigation of what quantities precisely 
are observable in the relativistic setting. Some observ- 
ables have been investigated previously, most notably the 
number density of tracers and the magnification. Here, 
wc have presented a unified relativistic analysis of "stan- 
dard rulers" , where a standard ruler simply means there 
is an underlying physical scale which we compare the 
observations to. This treatment applies to lensing mea- 
surements through galaxy ellipticities, sizes and fluxes, 
or through standard candles, to distortions of cosmolo- 
gial correlation functions, and to lensing of diffuse back- 
grounds. 

We show that in this framework, for ideal measure- 
ments, one can measure six degrees of freedom: a scalar 
corresponding to purely line-of-sight effects; a vector (on 
the sphere) which corresponds to mixed tranverse/line- 
of-sight effects; and a symmetric transverse tensor on the 
sphere which comprises the shear and magnification. We 
obtain general, gauge-invariant expressions for the six ob- 
servable degrees of freedom, valid on the full sky. These 
constitute the main result of the paper and are given in 
Eq. go]), dll]), (I51|), and ([64]). The vector component 
and the shear admit a decomposition into i?/_B-modes. 
The B-modes are free of all scalar contributions (includ- 
ing lensing as well as rcdshift-spacc distortions) at the 
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linear level, making them ideal probes to look for tensor 
perturbations. As an application of our results, we study 
the shear induced by tensor modes (gravitational waves) 
in 111. 

The logical next step is to construct estimators for 
these degrees of freedom, based on measurements of the 
density field of tracers (such as galaxies, the Lyman-a 
forest, 21cm emission, and so on). We will leave this for 
future work. 
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Appendix A: Spherical harmonic decomposition of spin-s functions 



Here we outline our notation and useful results on the spherical harmonic decomposition of tensors on the sphere. 
Throughout, latin indices ... denote components with respect to Euclidean coordinates, and are raised and lowered 
with Sij. We follow standard convention, see [32|. In particular, we do not include the Condon-Shortley phase in the 
spherical harmonics, so that (Y/m)* = i1-m- Explicitly, in our convention the spherical harmonics are given by 



(Al) 



where is a phase factor defined as 



1 , m > 
(-1)'" , m < 0. 



(A2) 



We can define spin±l unit basis vectors on the unit sphere 



^ I cos 9 cos (f> ± i sin ( 
m+ = —= {e\ T « e^) = —= cos sin (/) =F i cos ( 
v2 \/2 \ _sin6i 



(A3) 



where eg, are assumed orthonormal. m± transform as spinal fields (see §|T]). We have 



m!|_TO±i = 0, mt^m^i = 1 
•m^fii — 



(A4) 



Next, we define operators that raise and lower spin the spin s of a function (or tensor component) s/(0, </>) through 

lallai 



8J 



sin'* B 



i d 



d_ 

'do ^ sh^ 
d I a 



de sin 61 1 



sin-"0,J 
sin" 61,/. 



(A5) 
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A straightforward calculation using partial integration shows that 

J{d,f)s+ig dn^ Ji- sin^ 9) {[dg + t sin-i 90^] sm~' 9 ,/) ,+15 sin 9 d9 dcj) 
= J sfsin-"-^ 9{[de+i sin~^ 9d^] sin'^+i 9 s+ig) sin 9 d9 d(f> 
= J ,f{~B*s+ig)dn, (A6) 

in other words —8* is the adjoint operator of 9 with respect to the standard measure on the sphere. In many cases, 



which applied twice straightforwardly yields 



2.f{P; (t>) 



- e™^/(Ai), w 


2)(l+s)/2 


' d 


2)(l-s)/2 


■ d 


/ d 


m 






{--- 


m 







(A7) 



We can use the spin-raising/lowering operators to define spin-weighted spherical harmonics through 



{i + \s\y. \ (-1)^81^1 r,„„ s<o. 



(A8) 



(A9) 



Note again that our spherical harmonics are defined such that [Vim]* ~ Yi^^. Eq. (IA5p together with Ec^. (|A9p then 
yield [sYi„,]* = {—ly ^sYi~m- For s ~ ±2, this is equivalent to the definition used in |32l. ISq : 



±2Ylm 



G + 2) 



(AlO) 



The spin-weighted spherical harmonics defined in Eq. (|A9P form an orthornomal basis for spin-s functions. Using 
that 8^ = —8*, the orthonormality implies 



(/- \s\)\ 

Returning to the decomposition of a general spin±s field (with s > 0), we can express the components as 

±sA{-n) = ^at,±sYim{-^), 

Im 



(All) 



(A12) 



Acting with the spin-lowering operator on +sj4 and vice versa then yields 
8%A(n)= 



'■Im 1 



3W,™(n) = ^a,i(-l) 



S^_.A(n)=^ai(-l) 



Irn 



fcill8^e 
(i + sy 



'Im 1 





-sy 


{1- 


sy 




-sy 




sy 



Ylrnin) 

Yhn{n). 



(A13) 



We thus have 



J ±sA{h) [±,Yirain)]* dn 

l^^-^'^^\{-iy I [8^ +,A(n)] Y;^{h)dn 



ii + \s\y' 



' ii-\^\y- 
{i + \s\y 



[d^ ^sAinyY;jh)dn. 



(A14) 
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This shows that the coefficients have the desired property of being invariant under a rotation of the coordinate 
system around n. We can then define E- and i?-components through 



A _ AE . ■ AB 
AE _ 1 C^A , A*\ 



Under a parity transformation (n — > n' = — n), Y;,„ (— l)'yim- Hence, Eq. (jA9p yields sYi„i {—lY^sYir 
{—iy'^''[sYi_m]*, and E- and B- components transform under parity as 



af„^ ^ (-l)'+-'af_t 

af,^ ^ - (A16) 

Thus, the i?-component coefficients transform as expected of a spin-s quantity derived from a scalar perturbation; for 
example, a vector given by a gradient Bi = d±if, whose components ±iB transform as spin±l fields, is parity-odd 
(just like the electric field). On the other hand, the _B-component picks up an additional sign (parity-even, just like 
the magnetic field). 



1. Angular power spectra of spin-s functions 



The general procedure for obtaining the spherical harmonic coefficients and angular power spectrum for a spin-s 
quantity is as follows. The starting point is a relation between ^^(n, i) and the metric perturbations integrated 
over the (unperturbed) past light cone: 

,^(n, S) = r dx G[x, /i(xn, v)i (A17) 
Jo 

where rj = rjQ — /i(x, 77) stands for a single polarization state of any given metric perturbation, and the kernel G is 
a function of x ^ ^-nd its derivatives. Throughout, we suppress the polarization index, although all polarization 
states of course need to be summed in Eq. (|A17[) . The goal is to derive the contribution of a single Fourier mode of 
the metric perturbation, and to subsequently add up the contributions of all Fourier modes. For a scalar quantity 
(s = 0), such as a density or temperature, this calculation is straightforward since a scalar is invariant under a general 
rotation of the coordinate system. Thus, we can always align a given Fourier mode with the z-axis before summing 
up the contributions. For a general spin-s quantity gA, this is not possible. However, we can use the spin raising 
and lowering operators defined above to create a scalar quantity S'* (for s > 0) which allows us to easily sum up 
the contributions of different Fourier modes. The results of the previous section then immediately tell us how the 
resulting spherical harmonic coefficients of B'^ gA are related to those of gA. 
In detail, the calculation proceeds as follows: 

1. Evaluate the contribution sA{n,'k.) from a single plane-wave perturbation with a single circular polarization, 
with wavcvector k aligned with the z-axis. sA(n, k) is a function of n, usually written in terms of /x = ii • z and 
azimuthal angle 0: 

gA{h, k) r dx G(x, k, ^i);^(k)e-*e"^ (A18) 
Jo 

where /i(k) is the Fourier amplitude of the mode at some reference epoch, x = kx and r is an integer. In 
particular, r = if /i is a scalar metric perturbation, r — ±1 for a vector perturbation, and r = ±2 for a tensor 
perturbation, depending on the polarization state. G is an ordinary function obtained from G[x, h] by replacing 
di with iki, and pulling out h(k). Note that G thus contains the transfer function of the metric perturbation, 
relating h at the reference epoch to h at conformal time rj = rjo — x- 

2. Apply spin-raising or lowering operators to obtain a scalar quantity, 3* sA(fi,k) if s > 0, S'**' sA{h, ICj if s < 0. 
This quantity is a scalar on the sphere. By virtue of the exponential e*^'^, we can turn derivatives with respect 
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to fi into powers of ix, and powers of /i into derivatives with respect to ix. We can then write 

8^ sA{n, k) = rdxY] W,ixmx)il - ;,2)kl/2g.r0g..M /,(k)^ (A19) 
Jo i 

where Qi{x) are derivative operators in x, Wi{x) are coefficient functions, and we have pulled out a factor of 
(1 — for later convenience. Note that since the Qi{x) are constructed out of powers of ik = ix/x and 

d/d(ix), the terms involving even powers of x, dx are real, while those involving odd powers are imaginary. 
Hence, Q*{x) = Qi{-x). 

3. Since the angular dependence is now entirely in the factor (1 — /i^)l''l/^e*'''^e'^'', we can straightforwardly expand 
this scalar quantity in terms of the standard spherical harmonics Yim following Eq. (jAl4[) . The following relation 
which we prove in App. IA2l is useful: 



d^l y,:„(l - ^2)kl/2g..0g... ^ V47r(2/ + l)^/|^S|,'- (A20) 
where r is an integer. With this, we obtain 



at{k) = J^.VM2? + l)W j| jj + |^|j; (-ir dxY^W^ixy Mk). (A21) 



{l + \s\)\{l-\r\)\'- ' 7„ ' xM 

4. Following Eq. (|A16p . we can now separate the E- and i?-mode contributions 



AE 
'•Im 



(fc) ^ ^„,.VM2/ + l)^ ||^j^j|;|| + j^j|; (-l)- ^' £ dx J2 W^ixV Re Q,:(x)^ ^k) (A22) 
(fc) = ^,„.VM2/ + l)W ||~j^j|;|| + j^j|; (~l)- ^' £ dxj] M^.(x)*'- Img,(x)^ hik). (A23) 



These constitute the multipole coefficients of E- and i?-modes of sA. 

5. The angular power spectra are straightforwardly obtained by taking the expectation value of quadratic combi- 
nations of a^j^(fc), where X = E,B, summing over m, and integrating over {2'K)~^(Pk cPk' . Note that since we 
derived the multipole coefficients through the scalar quantity S'' i5^(n, k) which is invariant under rotations of 
the coordinate system, we can always align the Fourier mode with the z-axis, so that Eqs. (jA2ip - (jA23[) remain 
valid. 

d^k r d^k' 



rn— — l ' 

2 (^-|^l)!a+M)! 
7r(Z + |.|)!(/-|r|)!- 



\np J k'dkPh{k)\Fi''{k)\^ (A24) 

FP{k) = jyxY.W^{x)l^Q{kx)^ry (A25) 

Here, Ph{k) is the power spectrum of /i(k) at the chosen reference epoch, Np denotes the number of polarization 
states, and we have assumed that the different polarization states have independent phases and equal power 
spectra. 

Eq. (jA24p is a general expression for the i?/i3-mode angular power spectra of a spin— s observable induced by a spin— r 
metric perturbation, which is straightforward to evaluate once an expression of the form Eq. (jA19p is given. Note 
that for r = s = and Np = 1, we recover the usual result for scalar obscrvables induced by scalar perturbations. In 
the following, and in the related papers , we will apply this result for s = 0, ±1 and ±2, as well as r = 0, ±2. 
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2. Proof of Eq. (|A20)) 

The useful relation Eq. (|A20p is easily proven by induction over r. First, the case r = 0, 



dn YC^e'-^^ = V47r(2/ + l)i[n (x)(5„o, (A26) 
follows immediately from our definition of the spherical harmonics, Eq. (|A1[) . and the partial wave expansion 



^(2Z + l)i'jz(a;)-P/(cos0). 



(A27) 



Further, the definition of spherical harmonics, Eq. (jAip . yields 



where 

Prix) = (-i)'"(i - x2)W2_p,(^) > 0) 

is the associated Legendre polynomial. Comparing with Eq. (jA20[) . the conjecture to prove is thus 

JrL . _ il + \r\y- M+iji(^) 
' ^ ' - {l-\r\)l' 



(A28) 



(A29) 



(A30) 



Wc now proceed the proof by induction, assuming that Eq. ()A30P holds for some r > (without loss of generality). 
Using one partial integration on I['^^{x) we obtain 



= i /■ dfi{l - m')('-+i)/2(-1)('-+i)(1 - 



-1 
1 



[(1 - Ai^)'^+ie*^^] 



dM-ir(i-A^'r 



[-2fi{r + 1) + {I - fi^)ix] e- 



[2{r + l)d. + x{l + dl)\ il[ix) 



(A31) 



where in the third line we have converted powers of i/i to derivatives dx- Now we use that, by assumption, Eq. (|A30[) 
holds for r, which leads to 



rr+l/ \ _ + ^V- -r+l+l 



^r+i+i \^2{r + l)dx + x{l + dl)] 



(A32) 



Straightforward algebra, together with the differential equation satisfied by spherical Bessel functions, = —2i'Jx- 
[1(1 + \)lx^ - leads to 



([+I±l)l,,.+i+iJiM 
(l-r-\)\ x'^+i' 



which proves the conjecture Eq. (|A30p . 



(A33) 

□ 



Appendix B: Geodesic equation and displacements 



In this appendix, we outline the derivation of the displacements Eqs. (|34[) - (P7|) in the general gauge given by Eq. ^ 
and Eq. (|19p . This is a generalization of the derivation in ^l| (who considered synchronous-comoving gauge), and a 
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special case of the treatment in [331. Choosing the (zeroth order) comoving distance as afhne parameter, the photon 
momentum can be written as 

^ = {-l + Si^,n' + Se'). (Bl) 
dx 

The geodesic equation then becomes 

d dx^" dx" dx^ _ 

dxdx^'-^dx dx ^""'^ 

The zero-th order parts just yield dx^ / dx =const. We now turn to the first order part. The temporal component 
gives 

d ^ n dx°' dxl^ 



±5v + A'-2A..,7i}+(]-^' 
dx V2 



5v + A' - 2A,in' + ( -h'^^ + ) fi'h^ = 



|-(5:.-2A) = A'-i^-a||S||. (B3) 



The spatial components yield 



d , , dx"" dx^ ^ 

-^Se^ + A-^ B^' 2\h')n^ ~ 2\{B,, - B,^,)W + \ (h^ + h\^^ h^^ f^f^ = 

A (5e' + + h\jh^) = -A'' + ^^B\\ - Si, + i/^ii ' - ^V'^h.kfi''. (B4) 

Next, we need to obtain the initial conditions at the observer for the quantities 5e*. For this, we consider an 
orthonormal tetrad (e°)^i, defined through 

5'"^(e'^)p(e^). = (B5) 
To zeroth order, we set (e'')^ = a^Ty^^. Using Eq. ([2]), this yields at first order 

-1 = .g^''(e°)^(e"), = a-\-l + 2A)[{e%f 

= .9^''(e°)^(e^). = [-(e°)o(e')o - S-''(e'),(e")o + (e°),(e*F] . (B6) 

The first line yields (e'')o = a(— 1 — A), while the second line yields (e'^'); = a ((5'^; + h\/2y With this, the third line 
becomes 

= a(e')o+a^B' + a(e°),. (B7) 

Further, we require that the spatial hypersurfaces defined by (e*)^ be orthogonal to the four-velocity of comoving 
observers (Eq. ([T9)) ): 

= (e')^ u^^ = (e')o a-\l - A) + v\ (B8) 
Hence, we have (e')o — — au*, and Eq. (|B7p yields (e°)i — a{vi ~ Bi). We thus have 

(e°)^ =a(-l-A, v,-B,) 



(B9) 



Comparing the first line with Eq. (fT5|) . we indeed see that (e")^ — w^, as desired. 
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We now require that the normahzed photon momentum at the observer, pi^ = (—1 + SvotTi^ + ^e* ), measured with 
respect to this tetrad has the components (1, n*) (note that for our choice of afhne parameter, the photon momentum 
is past-directed). Setting Oo = 1 without loss of generahty, we obtain 



1 = (e°)^ p^' ^l + Ao-Siyo + {v.,o - B,o)i 

1 

2' 



fe = (e^)^p^ - + + + -(h^^of^. 



(BIO) 



Note that the first condition imphcs that a locally emitted photon is observed with unperturbed rcdshift (see 
Eq. (|BT6l) V Hence, 



Eva = Ao + V\\o - B\\o 



5el 



(Bll) 



Note the aberration term — in (JeJ,. We can now integrate the geodesic equations given these initial conditions: 



5v{x) = 2A{x) + / dx! 



Ao + v»,-B», + 2A{x)+ / rfx' 



const. 



The spatial component yields 

Se^x) = \{h\^o - + - B\x) - h'jix)n' + j dx' 
Integrating again yields the temporal and spatial displacements: 



~A' + h^Bj,, + -h^^''7Vn'' 



(B12) 



(B13) 



Sx'ix) 



Sx\x) 



dx dv{x) 

ri 
Jo 



dx 



2A{x) + {x^x)\A' --h^-d^^B^ 



dx^e\x) 



dx 



B' - W^n? + {x~x)\ ~A'' + fiPBj,, + -h^^f^fi^ 



(B14) 



(B15) 



Next, we need to evaluate the scale factor and affine parameter at emission, by requiring that the observed photon 
frequency match the redshift z. Since the photon momentum is given by Eq. (|Bip . we have 



l + z = 



{a-^u^dx''/dx)e _ a-^{x°){l + A - Sv + -t;|| - S||)e 
{a-^Uf,dxf^/dx)o 



{l + A-Sv + v^^ - B\\)o 



(B16) 



where a subsript e denotes the emission point, and we have decomposed Bi ~ B^^fii + B±i. The initial conditions for 
Sv imply that the denominator is 1, and we obtain 



aix°) 



1 + A- 5v + v» - B» = 1 + A In a, 



(B17) 



where all quantities on the right-hand side are evaluated at emission, and we have defined the perturbation to the 
logarithm of the scale factor as emission, A In a. Explicitly, 



A In a = A - + U|| - B\\ 
= Ao- A + v\\ - v\\o 



dx 



-A' + -h[ + B\ 



(B18) 
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Since a = 1/(1 + z), and a{x^) = 1/(1 + 2;) where z is the redshift one would observe for the same source in an 
unperturbed Universe, we can write this as 



l + z= (l + z)(l + Ahia), 



(B19) 



The perturbation to the conformal time at emission has two contributions, from the temporal displacement 5x^ and 
from the perturbation to the afhne parameter from its zeroth-order value, x'- 



Eq. ((BT9| thus yields 



so that 



, ns 91na r n „ n 

a(a;") = a + [5x" - 5x\ ■ 



l = l-hH{i) [5x° - (5x] + A In a. 



5x (^x" — ^rrrrr A In a. 



H{2) 



We can now assemble the total line-of-sight deflection: 



l + z 



Axil = Sx^fii + Sx = Sx\\ + Jx" ^Alno 



dx 

Jo 

For the transverse deflection, we obtain 
Aa;^ = dx^ — n'Sx\\ 

\v''{h,k)o + Slo - «lo] X+ r dx 



1 + ^A1 



(B20) 



(B21) 



(B22) 



(B23) 



(B24) 



which can be further manipulated to yield Eq. p6p . 



Appendix C: Test cases for the shear 

In this appendix, we consider test cases in order to validate the expression for the shear, Eq. (|65p in § IVII Bl For 
a larger set of test cases applied to scalar quantities such as the observed galaxy density and the magnification, see 
Appendix C in [llj- 

The first case is a metric perturbation corresponding to a pure gauge mode, i.e 

h,j (x, ri) = + Byfea;'=, (CI) 

where Aij and i?yfc are constant and symmetric in i and j. Such a metric perturbation can be obtained through a 
coordinate transform .t* .t* = — AijX^ /2 — BijkX^ x^ / A. Choosing the observer to be at the origin, we have 



L 



h±s =A± + B± kx" 



^ dx 1 



X ^ " 2 

(C2) 



where B±k = rn^j-mi'-^Bijk.. Since didjhki = and dihki = Bku = const, Eq. (|^5|) then yields 



±27(n) = -\a± - \^A± - ^B± ki'' + A± + ^B±k£'' = 0. (C3) 



19 



Thus, a gauge mode does not contribute to the shear. 

Further possible test cases are a perturbed expansion history and spatial curvature. Both cases correspond to 
isotropic Universes, and the observed shear should thus be zero. The former case is described by a metric perturbation 
of the form hij = A{T])dij. Since this implies that dihjk = and h± = 0, Eq. (j65p implies ±27 = 0. Spatial 
curvature is decribed by a metric perturbation hij = —K/Axkx'^Sij. In this case, h± = 0, dih^i = —K/2 XiSki, and 
^^^Jhkl=-K/2S^JSkl. Thus, 

K 

[ml^m'^didjhM)n^h^ = - — {m'^m^if = 

n'^m^m^dihki = ~^{m^Xt)hkm'^ = 0. (C4) 

This implies that ±2"/ = for spatial curvature as well. There is however one test case where the shear is non-trivial, 
which we will consider next. 



1. Bianchi I cosmology 



A Bianchi I cosmology is an anisotropically expanding Universe. Following |ll| . we choose the 3-axis to be unper- 
turbed, while the scale factors along the 1- and 2-axcs arc perturbed in the following way: 



aiiv) = a.{r])[l + si(r/) - 53(7;)] 
a2{il) = a(?7)[l + S2(r/) - 83(1])] 



(C5) 



where si{ri) + 32(77) -I- 33(7]) = 0, and 5^(770) = 0- Relaxing either of these conditions leads to cases we have studied 
above (perturbed expansion history and pure gauge mode). The non-zero components of hij are then given by 



/ill =2[.si(7y)-S3(7,)] 

/122 = 2[S2{7]) - 83(7])]. 



(C6) 



Let us consider two lines of sight close to the unperturbed 3-axis. Specifically, we consider photon 4-momcnta at 
the observer given by p'j}^ = (—1, — ^jO, —1) and plf-* = (—1,0, — ';,0, —1), where <r is the infinitesimal angle with the 
3-axis. By following back these geodesies, one can straightforwardly derive the angular diameter distances along the 
3-axis, for an object extended along the 1- and 2-axes 



DA^physAv) = a('y)(% - v) 



1 + Saiv) - S3iv) + 



(C7) 



where a = 1, 2, and rj is the conformal time of emission. The observed ellipticity of galaxies, designed as estimator for 
shear, can be written as 



1/11 



ID 



A.phys.l 



-D 



-2 

A,phys,2 



£2 



2 hi + I22 

/12 ^ 

hi 



2 



D 



-2 

A,phys,2 



0, 



'22 



where hj are the quadrupole moments of the galaxy's light distribution, which scale as Dj^^^^^. 



(C8) 

Other definitions 



are possible, however all of them agree at linear order. 62 vanishes, since a ray with = (— 1, — 1;/-\/2, S^/v^, — 1), 
propagating at -1-45° azimuthal angle to the 1-axis, yields the same angular diameter distance as a ray with p^^ = 
(— 1,— <;/-\/2, — <;/-\/2, — 1), propagating at —45° angle. Assuming that the galaxies are on average round {{ci) — 0), 
i.e. that they are not directly influenced by the anisotropic expansion, we obtain 



1 



1 



4 

2si + 2s3 + 2s2 - 2s3 / drj' [-si{t]') + 33(7/)] 

Vo - ?? J, J 
1 rvo 

„ .^2(7;) - Sl(7y)] + / d,/ [51(7?') - ,S2(7/)] . 

2 7/0 - V Jr, 



(C9) 



20 



The extra factor of 1/2 is due to the sum of moments in the denominator in the definition of ei. In order to compare 
with Eq. (|55|). we use 

h± = ^{hii - h22 ± 2ihi2) 

==si(r;)- 53(77). (CIO) 

We thus have h±o = 0, and Eq. (|65|) yields for the shear along the 3-axis 

1 dx 
71 ±^72= - 7;isi{v) - S2iv)) + / — [siivix)) - S2ivix))] 
^ Ja X 

= Mv)-si{7l)) + / dr/[si(7y')-S2(r]')], (Cll) 

where we have used that ±27 = 7i ± «72 for the coordinates chosen here (see Eq. ([50])'). Since this expression is real, 
72 = 0, and the result is equal to 71, which moreover agrees with the correct physical result for ei, Eq. ()C9|) . 



2. Shear from scalar perturbations 

In the derivation leading to Eq. (j65p . we have not made any assumptions about metric perturbations except that 
Sgoo = = Sgoi (synchronous-comoving gauge). As a cross-check, we now consider the case of scalar perturbations, 
where we can write 

h^j = 2D6,j + 2 (d,dj - ^-^yV^^ E (C12) 
in terms of the 3-scalar perturbations D and E (see [HI). We thus have h± — 2m^m^£'_ij and obtain 
±27(fi) = - m^m^ {E,,,{o) + E,,,{s)) - J^^xi^ix- x)jrn^rn^ (^D^,, - ^V^E^,, + c>f[S,,^ 

+ 2 ^1 - 2^^ m^m^d\\E^,^ - ^m'^m^E^,^^ . (C13) 

We now consider the contribution of a single plane wave along the z-axis, 

D{x, 1]) = D{k, T])e'^-^^ = ^(k, 77)e"^, (CM) 
and similarly for E. We can then replace 

m'^m^^^^J ^ -^(1 - M^)fc^- (C15) 



This yields 



±27(k,n)= +i(l-Ai')P(i?(k,77o)e"'' + E{k,f^y-^ 



x=0 



+ I dx{{x-x)jl{i-f^^)k'-(^D{k,rj) + h'-E{k,rj)-fi'eE{k,,^) 

1 - 2^] 7(1 - i^L^)fik^E{k,r]) - k^-{l - fi^)E{k,ij) le"^, (C16) 
Xj X ) 

where fj = tIq — X, and x = kx- Note that a scalar perturbation produces equal amplitudes of ±27^ there is no preferred 
handedness for scalar modes. Correspondingly, since this expression is oc e*™''^ with tti = 0, the spin-lowering and 
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spin-raising actions [Eq. (|A8|) ] become equivalent. We have 
e227 = 9'-27=|^[(l-M')27(ii,k)] 



(C17) 



x—kx 



Using that 



1-2^^ -QUx)k'E{\<i,i^)~-Qf{x)k'E{\^,r^) U^e"" 
Xj X X 



^ [(1 - ^^)^e«^] =^ [(1 + dlfe^^'^] = -(1 + a^)^ [x^e^^-] , 

1^ [(1 - = = ^9,(1 + [^^e^^-] , 

1^ [(1 - ^2)V^e"-] = = 9^(1 + , (C18) 

we obtain the operators Qf{x) for the scalar case: 

Qfix) = -{1+dir 

Qlix)=il + dlfdl 

Qi{x)=x{l + dl)d,. (C19) 

Note that all these operators are real. Hence, following the general derivation in App. lA 11 there are no parity-odd 
terms in the scalar contributions to the shear, and thus no i?-modes as expected. 



Appendix D: Vector 

In this section, we follow the procedure described in App. lA II in order to derive the multipole moments of the 
vector component ±iB induced by scalar and tensor perturbations. For scalar perturbations, we work in conformal- 
Newtonian gauge; there is no gauge ambiguity for tensor modes. Our goal is to show that scalar perturbations do not 
contribute to afj^, while tensor perturbations do contribute. 



1. Scalar modes 

We begin with Eq. 

(±i6)cN ^-v± + (^-d±^ + a±[i;|| - + dx|a±($' - vj/') j , (Dl) 

and specialize to the case of a single plane wave along the z-axis. Further, we write Vi = Vi. Noting that 

d±iV\\ = d±i(n^ diV) ~ h^d^idjV + —di_iV 
II ' X 

^L^v\\, = \{d^,V)o (D2) 
X 
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this yields: 



(±iS)cN(n,k) = -^fc±F(k,77)e"'' + 



rdx^*fc±[$'(k,7;)-*'(k,r,)]e-^ 
^0 X 



72 



X J X 



x=0 



''*^ia;[$'(k,7,)-*'(k,r/)]e'^^ 
X 



where as before x = kx, x = kx, and fj = rjo — X- Further, we have used that fcy = n*fci = /zfc, and k± = 
ml^ki = — a/1 — /i^ k/y/2. We now apply the spin-lowering operator Eq. (jA7P to iB. Note that since ±ii3(fi, k) is 
(/>- independent, this is identical to the expression for 9 



-B(n,k)=5^ VW^i6(n,k) 



zx[$'(k,7;)-*'(k,77)]e^ 

^0 X 



v2x 



V2i/(S) I X 



1 r 

X 



Qbs2{x) - -t^Qbsi{x) 



dx, 



QBSi(x)e"'' y(k,r?o)-/ ^QBsi(x)[$'(k,r?)-vI/'(k,77)]e'-^ 
^=0 ./o X 



(D3) 



where wc have defined 



X — 2xdx — X 



Qbs2{x) = I + xd.,. + idl + xdl. 



(D4) 



This expression is in the desired form, Eq. (jA19p . and we sec that all coefficients are real. Thus, scalar perturbations 
only contribute to the (polar) £^-coniponent of S^. 



2. Tensor modes 



We begin with Eq. (|45 



(±lS)s 



1 + 5 



1 + 5 



2H{~z)J^ " 2i/(i)7o ^x 



dXTd±h\\ = I dx% [ {d±hii)h'n' + hi'^im^^n^ 



Wc then decompose /ly into Fourier modes of two polarization states (see [2l|, [l^ for details) , 

/i.,(k,r?) =e+(k)/i+(k,r;) +ef,(k)/i"(k,r;). 



(D5) 



(D6) 



where e^j(k), s = +, x, are transverse (with respect to k) and traceless polarization tensors normalized through 
e^jC'* = 26'"^ . We assume both polarizations to be independent and to have equal power spectra: 



(/i,(k, ri)hs. (k', 7/)) = (2^)35p(k - k')(5,s.' -Trik, vWK v')PTo{k), 



(D7) 
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where Txik, rf) is the tensor transfer function, and the primordial tensor power spectrum is denoted as PTo{k). Further, 
we can define hehcity±2 polarization tensors and Fourier amplitudes through 



(D8) 



Note that Ph^^ik) — PTQ{k)/8. As before, we begin by evaluating the contribution of a single plane wave, assuming 
that k = fcz. We have 



1 



e5,(k,n) = ef^-(k)m*^(n)m^^(n) = -(1 TPm) e 



where fi = cos0. Restricting to the polarization p = +1 first, we have 



tiS(k, n,+l) 



1 + S 



X 



V2 



(1 _ ^,2^,k + ^^Y^i^L T 1) h[{k, r,) 



X 



(D9) 



(DIO) 



We now apply the spin-lowering operator Eq. (jA7p with m = 2 to obtain 

2 



giS(k,n,+l) = 



where 



1-^2 

2 

1-^2 



Vl-fihBik, n,+l) 



l + z 



2H{S)Jo X 



dx\ (1 - 



2\2 



V2 



-ix - V2(l - ^?){\ - ij) i /i'i(k, 77)e'2*e"'' 



l + z 



dx r 1 



2i/(5)7o X IV2 



(x-) + y2QBT2(a;) \ (1 - ^^2)g^20g^xA. U'i(k,r,) 



Qbti{x) = + Axd^ + ar^a^ + 2ia; 
Qbt2{x) = 3 + xdx - ix. 



(Dll) 



(D12) 



For the other polarization state, we obtain the same result with /i — > — x — > — x, </> — > —(f). Eq. (|D11[) is in the 
desired form, Eq. (jA19[) . We see that the operators have both real and imaginary parts, signaling that tensor modes 
contribute to both the polar ("i?-mode", through ^cQBTi) and axial vector ("_B-mode", through ImQsTi)- 



Appendix E: Components in terms of matter density contrast 

This section gives useful expressions for C, ±iB, and the shear in terms of the familiar matter density contrast 
(5^ in synchronous-comoving gauge. For convenience we write the velocity Vi in terms of a scalar velocity potential 
V, Vi = V^i, and relate V, $, ^E* to the density contrast (5^ in synchronous-comoving gauge through (see p38l. l33|) 

F(k,r;) =ai//fc-2i?(a(,7))<5^(k,,yo) 
$(k, 77) - *(k, ry) = (fc, ry)C(k, 7?o) 
$(k, 77) + vl/(k, 7?) = .g(fc, 77)Z?$_ (fc, ?7)C(k, ?7o) 

^$(k,r;)=i[.9 + l]^*-C(k,?7o) 

*(k,77) = i[5-l]i?*_C(k,77o) 
($ - vl/)'(k, r;) = Asw(fc, ry)C(k, m), (El) 
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where / = dlnD/d\na, D{a) is the matter growth factor (normahzed to unity at a = 1) and we have defined general 
coefficient functions -Disw to allow for non-standard cosmologies. In a ACDM cosmology (or more generally 

for a smooth dark energy component) , we have 

5(k,77)=0. (E2) 



Here, a subscript denotes that the quantity is defined at the present epoch — rjo. We will denote the power 
spectrum of (5^ at z = as Pm{k), as well as x = kx, x = kx as before. 



1. Longitudinal scalar 



We begin with Eq. (gj), 



(C)cN = - fe^Alna-^-wii + + - + 



(5) 



(E3) 



where 



1 + z \ 9 In ro 



H{z) J 9 In a 



Alna = *o - * + I'll - I'll 



dx W - ^-1 . 

JO 



(E4) 



Note that for a non-evolving ruler in ACDM, hz ~ 3f2m(z)/2. For a single plane-wave perturbation, this yields 



C(k, n)^hA[ *(k, 77) - ^^Jix^^^^ I e 



x=0 Jo 



dx [$' - *'] e"'^ 



X 



1 + z 

H{~z)x 



VCk fj) 



1 , , aHfD - 

2 A:^x 



*_ rT-^Qcsi[x) 

k^X 



dxDiswe 



+ 5-(k,7yo)j- 



1. A aHfD 



HCz)x 



^{g - l)D^_Qcsi{S:) + ^^i^Qcs2{x) - [aH fD)'Qcsi{3:) 



k'^X 



+ I {g'D^- + (.9 + i)Asw) 



k,r} 



(E5) 
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where Qcsi{x) = xd^, and QcS2{x) = x^d^. This is clearly in the form Eq. (jA19|) . with r = s = 0, and we can thus 
immediately apply Eq. (|A24p : 



Cc(0 



edkPra{k)\Ff{k)\'' 



Fi:{k) = b.Fr'"'{k) + { 



Fi^^'^'ik) 



-{g-l)D.^_+-^Qcsi{x) 



H{i)x 



i(5 - l)D^_Qcsi{x) + ^^^Qc52(i) - {aHfD)'Qcsi{i) 



X 



{g'D^_ +{g + 1) Asw) 



3i{x) 



k,ri 



aHfD ^ I, , \ , , 

— - 77(5 - 1)-D<i>_ ji(x)+ dxDiswji{x). 



(E6) 



2. Vector 



As derived in App. [Dl 



V2x 



V2H{i) [ X 
1 



?BS2(S) - ^QBSl(i) 



^2 



^=0 Jo X 



aHfD 



V2k'^ 



X 



)bsi{x) 



l + S 



V2H{S) \2x 
l + S 



— [g- 1]D^_ (fc, a)QBSi{i) + aHfD 



)bS2{x) - -^QBSlix 



V2H{i) 

Since ±iB is a spin-1 quantity, Eq. (IA24P with r = 0, s = 1 yields 



, ''?^Asw(fc,a(x))Qs5i(x)e"^ 

2:=0 Jo X 



(E7) 



2 1 



7rZ(/ + l) 



k'dkPrn{k)\FP^{k)\^ 



Frik) = 



zQbsi(x) 



l + z \ 



V^k^X 



V2H{z) \2x 



-[g- 1]D^_ [k, ~a)QBSi{i) + aHfD 



QbS2{x) - -^QBSlix) 



{HfD),=o^ QBSi{x)ji{ 



—Disw{k,a{x)) QBSi{x)ji{x) 
^=0 Jo X 



(E8) 



This is the power spectrum of the i?-mode (polar) component of ±B. The operators Qbsi as applied to spherical 
Bessel functions become 



QBSi{x)jiix) = + l)jiix) 
Qbs2{xMx)= 1(1 + 1)8, ^^'^''^ 



(E9) 



Note that this implies [QbS2{x) — x '^QBSi{x)]ji{x) = 1(1 + \){dxji{x)) / x. Note further that the observer term 
\QBSi{x)ji{x)]x=Q is only non-zero for the dipole / = 1, as expected. 
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3. Shear and magnification 

We begin with the shear. Eq. ((66| evaluated for a plane-wave perturbation yields 



Jo A ^ 



dxix- x)jk^D^_ {k, a(x))(l - /i')e'^^ 5^^(k, 7?o). 



In order to obtain a scalar quantity, we now apply Eq. (jASp with to = 0, yielding 
8227(n,k)=a^[(l-M^)27] 



(X - xnk^D^^ {k, a{x))dl [(1 - ^l^fe^-^] <5-(k, 
X 

rfx (X - x)^k'D^- {k, a{x))Q^s{x)e^''^ <5^^(k, ,70), 
X 



where we have defined 



Q^s{x) + + Sxda: + (12 + 2x^)dl + 8xdl + x^d^. 



(ElO) 



(Ell) 



(E12) 



This is again in the desired form Eq. (jA19p . and applying Eq. (|A24p with r = 0, s = 2 yields the angular power 
spectrum of the shear iJ-modes, 



2([-2)! 
^'^-^(7T2)! 



edkP^{k)\F^'^{k)\ 



1 



dx (X - x)^k'^D^_ {k, a{x))Q^s{x)ji{x). 
X 



(E13) 



The operator applied to spherical Bessel functions is 



Qis{x)ji{x) 



{l + 2)\ji{x) 
[1-2)1 x^ ■ 



(E14) 



The magnification contains a number of terms, which are evaluated in Appendix B of |24l |. Here we only repeat the 
final result: 

Cm{1) = ^ f k^dkPrn{k)\Fi^{k)\^ 



Fr{k)^2 



aHx 



Fr''%k)-{[g + l]D^_),]i{i) 



+ + I dx^^D^_ji{x)+2 I ^D^Jiix) 



XX 



dx 



X 



(E15) 



where Ff'^''°'{k) is defined in Eq. (|E6)) . 
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